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PREFACE 


This  report  was  prepared  by  Gary  Koh,  Physical  Scientist,  Geophysical  Sciences 
Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory.  Funding  for  this  research  was  provided  by  DA  Project  4A161 102AT24, 
Research  in  Snow,  Ice  and  Frozen  Ground,  Scientific  Area  A,  Properties  of  Cold 
Regions  Materials,  Work  Unit  004,  Properties  of  Cold  Regions  Hydrometeors.  The 
author  thanks  Harold  W.  O’Brien  and  Dr.  Richard  H.  Munis  of  CRREL  for  their 
suggestions  and  technical  reviews  of  this  report. 

This  report  describes  a  preliminary  investigation  of  the  extent  and  military  applic¬ 
ability  of  present  knowledge  of  ice  fog  as  an  electro-optical  obscurant.  Relationships 
are  developed  that  may  contribute  to  the  Electro-Optical  Systems  Atmospheric  Ef¬ 
fects  Library  developed  at  the  Atmospheric  Sciences  Laboratory,  as  well  as  provid¬ 
ing  initial  guidelines  (rules-of-thumb)  for  arctic  sensor  selection  and  utilization. 
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ICE  FOG  AS  AN  ELECTRO-OPTICAL 
OBSCURANT 

Gary  Koh 


INTRODUCTION 

Ice  fog  is  an  atmospheric  phenomenon  visually 
resembling  common  water  fog.  However,  ice  fog 
usually  occurs  only  when  the  ambient  temperature 
drops  below  -30  °C,  and  generally  increases  in  fre¬ 
quency,  intensity  and  persistence  with  decreasing 
temperature.  Under  a  microscope,  the  particles 
composing  ice  fog  are  seen  to  be  minute  ice  crys¬ 
tals,  approximately  2  to  20  /un  in  diameter. 

During  dense  ice  fog  the  visibility  may  be  re¬ 
duced  to  only  a  few  meters.  The  optically  degrad¬ 
ed  environment  associated  with  ice  fog  generally 
occurs  when  water  vapor  resulting  from  combus¬ 
tion  processes  condenses  on  particulate  matter  to 
form  small  droplets  that  rapidly  freeze  into  ice 
particles  when  introduced  into  a  cold,  stagnant  air 
mass  (Benson  1970).  In  inhabited  areas,  exhaust 
gases  from  home  heating  plants,  factories,  auto¬ 
mobiles  and  aircraft  are  notorious  contributors  of 
raw  materials  for  ice  fog  production. 

The  formation  of  ice  fog  with  resultant  low  visi¬ 
bility  can  seriously  restrict  the  effectiveness  of 
military  operations.  A  large-scale  military  opera¬ 
tion  conducted  in  arctic  and  subarctic  regions  with 
associated  vehicular  emissions,  gunfire  and  other 
combustion  activities  is  a  potential  source  of  the 
nucleating  material  and  water  vapor  required  for 
ice  fog  growth.  With  increasing  dependence  on 
electro-optical  weapon  systems,  a  model  predict¬ 
ing  the  performance  of  electro-optical  systems 
during  ice  fog  is  consequently  of  fundamental  im¬ 
portance  in  the  winter  battlefield  scenario.  This 
necessitates  a  detailed  knowledge  of  the  processes 
of  ice  fog  formation,  the  physical  and  optical 
properties  of  the  ice  fog,  and  an  understanding  of 
the  electro-optical  theory  involved. 


ICE  FOG 

Much  of  the  research  on  the  physical  and  opti¬ 
cal  properties  of  ice  fog  has  been  conducted  at  the 
Stanford  Research  Institute,  the  Geophysics  Insti¬ 
tute  of  the  University  of  Alaska,  and  the  U.S.  Ar¬ 
my  Cold  Regions  Research  and  Engineering  Lab¬ 
oratory.  The  properties  of  the  ice  fog  discussed  in 
this  report  were  obtained  from  the  extensive  infor¬ 
mation  provided  by  these  organizations.  The  pur¬ 
pose  of  this  review  was  to  extract  the  basic  con¬ 
cepts  relevant  to  the  study  of  ice  fog  as  a  cold  re¬ 
gions  obscurant. 

Effects  of  atmospheric  conditions 
on  ice  fog  formation 

The  effects  of  battlefield  activities  on  the  for¬ 
mation  of  ice  fog  strongly  depend  on  the  nature  of 
the  activities  and  the  ambient  conditions.  The  abil¬ 
ity  to  predict  the  impact  of  vehicular  operation  or 
gunfire  in  producing  ice  fog  is  of  importance  prior 
to  tactical  engagement.  Therefore,  the  environ¬ 
mental  conditions  necessary  for  the  formation  of 
ice  fog  as  a  result  of  combustion  activities  should 
be  identified. 

A  graph  (Appleman  1953)  of  the  temperature 
and  humidity  requirements  for  the  formation  of 
ice  fog  due  to  moisture  produced  by  the  burning 
of  hydrocarbon  fuels  is  reproduced  in  Figure  1. 
According  to  Appleman,  ambient  temperature 
and  relative  humidity  will  determine  whether  add¬ 
ed  moisture  and  heat  from  combustion  will  lower 
or  raise  the  humidity  of  the  affected  environment. 
The  critical  temperature  necessary  for  the  humidi¬ 
ty  to  increase  from  combustion  activities  as  a 
function  of  initial  humidity  at  1000  mb  is  repre¬ 
sented  by  curve  AB.  Thus  ambient  temperature 
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Figure  I.  Temperature  /humidity 
relationships  for  ice  fog  forma¬ 
tion. 


and  humidity  conditions  must  be  in  the  region  left 
of  the  curve  for  saturation  to  occur  from  added 
moisture  from  combustion.  As  an  example,  the  in¬ 
itial  humidity  must  exceed  50%  for  ice  fog  to  oc¬ 
cur  from  combustion  activities  at  -35  °C,  provided 
that  enough  moisture  to  produce  saturation  is  re¬ 
leased. 

Warming-up  and  subsequent  operation  of  air¬ 
craft  or  a  large  number  of  vehicles  under  the  ap¬ 
propriate  conditions  will  produce  ice  fog.  The  ex¬ 
tent  of  the  ice  fog  will  depend  on  the  amount  of 
moisture  present.  The  temperature  and  humidity 
regime  where  the  condensation  of  droplets  can  oc¬ 
cur  may  be  either  ice  saturated  or  ice  subsaturated 
(Fig.  1).  If  the  vehicles  are  operating  in  an  ice  sub- 
saturated  environment,  the  dissipation  of  ice  fog  is 
likely,  and  one  can  expect  localized  plumes  of  ice 
particles  around  the  vehicles  emitting  the  exhaust. 
However,  the  introduction  of  moisture  in  an  area 
near  or  at  ice  saturation  can  produce  ice  fog  that 
may  linger  and  whose  duration  will  again  depend 
upon  the  amount  of  moisture  introduced  by  the 
battlefield  activities  and  the  rate  of  ice  fog  settling. 

Ice  fog  is  associated  not  only  with  vehicle  opera¬ 
tion  but  also  with  military  activities.  Kumai  and 
O’Brien  (1964)  observed  the  effects  of  155-mm 
howitzer  activities  in  producing  ice  fog  and  reduc¬ 
ing  target  visibility.  The  physical  characteristics  of 
the  gunfire-generated  ice  fog  were  similar  to  those 
of  ice  fog  generated  from  other  sources.  Camp  ac¬ 
tivities  of  troops  in  a  cold  environment  are  an  ad¬ 
ditional  source  of  ice  fog.  The  sources  of  ice  fog 
are  important  to  identify  because  the  localized  ice 
fog  may  be  visible  for  many  miles,  with  the  poten¬ 
tial  of  indicating  the  troops’  presence. 


The  potential  exists  not  only  for  predicting  the 
formation  of  ice  fog,  but  for  determining  its  phys¬ 
ical  properties,  given  the  nature  of  the  water  vapor 
source  and  the  ambient  conditions.  The  general 
theory  of  atmospheric  ice  crystal  formation  found 
in  the  literature  of  cloud  physics  can  be  applied  to 
the  growth  of  ice  fog  if  the  different  environ¬ 
mental  conditions  associated  with  the  two  phe¬ 
nomena  are  taken  into  account.  A  theoretical 
study  on  the  formation  and  growth  of  ice  fog  has 
been  conducted  by  Huffman  (1968).  By  applying 
the  laws  of  thermodynamics  and  statistical  me¬ 
chanics,  Huffman  obtained  the  numerical  values 
for  the  size  distribution  of  ice  fog  in  general  agree¬ 
ment  with  his  experimental  results.  A  numerical 
model  by  Nelson  (1972)  for  computing  the  size 
and  concentation  of  ice  fog  particles  produced  by 
automobile  exhaust  again  showed  good  agreement 
between  the  experimental  and  computed  results. 
These  models  indicate  that  it  may  be  possible  to 
utilize  a  similar  computer  program  to  numerically 
characterize  the  ice  fog  resulting  from  battlefield 
activities. 

An  application  of  these  numerical  models  is  the 
determination  of  ways  to  reduce  or  enhance  the 
transmission  properties  of  the  ice  fog  generated  by 
vehicles.  Nelson  was  able  to  affect  the  concentra¬ 
tion  and  the  size  of  the  ice  fog  particles  by  varying 
some  operating  parameters  of  an  automobile  as 
inputs  to  his  model.  The  size  of  the  ice  fog  parti¬ 
cles  could  be  increased  by  increasing  the  exhaust 
pipe  diameter  or  the  exit  temperature  of  the  ex¬ 
haust  gas.  Similar  results  can  be  achieved  by  de¬ 
creasing  the  exit  velocity  of  the  exhaust  gas  or  the 
concentration  of  soluble  and  insoluble  particles  in 
the  exhaust.  Given  the  dependence  of  the  extinc¬ 
tion  properties  of  ice  fog  on  the  particle  size/ wave¬ 
length  ratio  and  the  size-dependent  settling  rate  af¬ 
fecting  the  duration  of  the  ice  fog,  further  investi¬ 
gation  into  numerical  modeling  may  be  of  interest. 

Size  distribution  of  ice  fog  particles 

Knowledge  of  particle  size  distribution  is  re¬ 
quired  to  calculate  the  attenuation  of  electromag¬ 
netic  radiation  transmitted  through  ice  fog.  For 
the  parametric  modeling  of  ice  fog  attenuation, 
the  substitution  of  experimental  size  distribution 
data  into  an  analytical  expression  is  desirable.  One 
analytical  expression  capable  of  realistically  de¬ 
picting  the  diverse  shapes  of  the  ice  fog  distribu¬ 
tion  is  the  modified  gamma  distribution  used  by 
Deirmendjian  (1969)  to  model  various  atmos¬ 
pheric  aerosols. 

The  modified  gamma  distribution  is  a  three- 
parameter  model  expressed  as 


Figure  2.  Modified  gamma  distribution  for  ice 
fog. 


n(r)  =  ar“  e\p[(-a/y)(r/rc)']  (1) 

where  n(r )  is  a  continuous  and  integrable  function 
of  radius  r.  The  shape  of  the  distribution  is  gov¬ 
erned  by  a  (positive  integer),  y  (positive  constant) 
and  rc,  which  is  the  mode  radius  of  the  distribu¬ 
tion.  The  value  a  normalizes  the  total  number  of 
particles  in  a  unit  volume  and  does  not  affect  the 
shape  of  the  distribution  curve. 

The  temperature  dependence  of  ice-fog  crystal 
growth  has  been  observed  by  various  investi¬ 
gators.  Thuman  and  Robinson  (1954),  Kumai 
(1964),  Kumai  and  O’Brien  (1964),  and  Ohtake 
(1970)  employed  impactors  or  electrostatic  precip¬ 
itators  to  collect  and  analyze  ice  fog  crystals  and 
noted  the  increase  in  the  mean  particle  size  with  in¬ 
creasing  temperature.  The  composite  size  distribu¬ 
tions  of  ice  fog  occurring  at  temperatures  of 
-31.0°  to  -32.9 °C,  -35.0°  to  -36.9°C  and  -39.0° 
to  -41.0°C  reported  by  Ohtake  and  the  corre¬ 
sponding  modified  gamma  distribution  fit  are  il¬ 
lustrated  in  Figure  2.  The  modified  gamma  fit  was 
approximated  by  trial  and  error  because  of  the 
limited  data  set.  Although  the  fit  does  not  neces¬ 
sarily  represent  the  optimum  values,  the  versatility 
of  the  three-parameter  model  is  demonstrated. 

For  future  studies,  optical  measurement  tech¬ 
niques  can  be  used  to  obtain  larger  data  sets  of  ice 
fog  size  distributions.  This  will  necessitate  a  com¬ 
putational  technique  to  obtain  the  best  fit  of  the 
parameters  to  the  measured  data.  Optimization 
programing,  evaluating  the  deviation  of  the  model 
distribution  from  the  experimental  data,  is  one  po¬ 


tential  method.  This  can  be  accomplished  by  mini¬ 
mizing  a  function  of  the  form 

m 

F\WC)  =  ^  [/>.  -/i(r,)P  (2) 

/  =  1 

where  m  is  the  number  of  experimental  data  pairs 
( m ,  ri),  and  n(r,)  is  the  modeled  distribution. 

Shape  distribution  of  ice  fog 

An  electromagnetic  wave  incident  on  a  particle 
is  scattered  in  all  directions.  The  intensity  of  the 
scattered  radiation  measured  at  a  given  location  is 
the  resultant  scattered  radiation  from  the  various 
parts  of  the  scattering  object.  The  phase  differenc¬ 
es  of  the  scattered  radiation  become  important 
when  the  size  of  the  particle  is  comparable  to  or 
larger  than  the  wavelength  of  the  incident  radia¬ 
tion.  Comprehensive  understanding  of  the  scatter¬ 
ing  process  requires  not  only  this  size  factor  to  be 
considered  but  also  the  shape  factor,  since  the 
phase  relation  of  the  scattered  radiation,  and  thus 
the  angular  scattering  pattern,  are  also  dependent 
on  the  geometry  of  the  scattering  object. 

The  classification  of  ice  fog  particles  according 
to  shape  conducted  by  previous  investigators  indi¬ 
cates  that  three  general  types  of  ice  fog  are  pres¬ 
ent.  Thuman  and  Robinson  (1954)  encountered 
mostly  droxtals  (near-spherical  particles  with  rudi¬ 
mentary  crystal  faces)  during  ice  fog  occurring  at 
temperatures  below  -38  °C.  However,  with  in¬ 
creasing  temperature  they  observed  that  the  per- 
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Figure  3.  Shape  distribution  of  ice  fog. 


centage  of  hexagonal  plates  and  column-shaped 
crystals  increased  rapidly.  Investigation  of  ice  fog 
formed  at  temperatures  below  -37  °C  (Kumai 
1964,  Kumai  and  O’Brien  1965)  also  indicated  the 
prevalence  of  spherical  particles  with  some  plate¬ 
like  and  columnar  structures.  Ohtake’s  (1970)  ob¬ 
servations  were  consistent  with  those  of  previous 
investigators.  He  reported  that  hexagonal  plates 
and  columns  accounted  for  over  60%  of  the  ice 
fog  crystals  at  temperatures  above  -37  °C.  The 
temperature  dependence  of  the  shape  distribution 
of  ice  fog  (from  the  findings  of  Ohtake)  is  illus¬ 
trated  in  Figure  3. 


ELECTROMAGNETIC  THEORY 
RELEVANT  TO  OPERATIONS 
IN  ICE  FOG 

Electromagnetic  radiation  propagating  through 
an  ice  fog  medium  is  scattered  and  absorbed,  thus 
incurring  a  path  attenuation  by  a  factor  of  exp[-fl.i, 
(X)L).  is  the  forward-scatter  corrected  ex¬ 

tinction  coefficient  of  the  medium  for  wavelength 
X,  and  L  is  the  path  length  of  the  incident  beam 
traversing  the  ice  fog.  This  is  the  familiar  Beer- 
Lambert  relationship,  and  the  term  enclosed  in  the 
brackets  represents  the  effective  optical  thickness 
of  the  medium.  The  extinction  coefficient  for  a 
volume  of  air  containing  ice  fog  particles  is  deter¬ 
mined  by  integrating  the  extinction  properties  of 
the  individual  particles  of  the  volume  and  is  ex¬ 
pressed  as 


'max 

B.„  =  J ■nr2n{r)[Q,cu(m,x)  +  Q.b,(m,x)\dr  (3) 

'min 


where  xr2  =  geometrical  cross-sectional  area 
n{r)  -  size  distribution  function 
(7*ca  »  Qofci  =  efficiency  factors  for  scattering  and 
absorption 

m  =  complex  index  of  refraction  of  ice 
fog 

x  =  size  parameter,  2xrA. 

The  scattering  efficiency  factor  in  this  equation 
fails  to  take  into  account  directional  scattering, 
and  thus  neglects  the  forward-scattered  compo¬ 
nent  entering  the  field  of  view  of  the  detector  dur¬ 
ing  transmittance  measurements.  Therefore,  in 
order  to  compare  the  theoretical  and  experimental 
extinction  coefficients,  it  becomes  necessary  to 
modify  eq  3  to  the  following  form: 

rmox 

B.„  =  J  *r2n(r){Q,c.-Q;ca  +  Qab,]dr  (4) 

rmin 


where  QL  is  the  efficiency  factor  associated  with 
the  scattered  radiation  entering  the  field  of  view  of 
the  detector  and  is  determined  from  the  angular 
scattering  properties  of  ice  fog  and  the  geometry 
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Figure  4.  Mie  efficiency  factors 
for  scattering  and  absorption. 


of  the  transmission  system.  The  multiple  scatter¬ 
ing  effects  are  not  considered  in  this  report. 

Because  of  the  mathematical  difficulties  en¬ 
countered  in  solving  radiation  scattering  by  ran¬ 
domly  oriented  asymmetrical  particles,  modeling 
of  electromagnetic  radiation  through  atmospheric 
aerosols  generally  assumes  the  sphericity  of  the 
particles.  Rayleigh  scattering  and  geometrical  op¬ 
tics  are  the  approximate  solutions  available  for 
transmittance  calculations  through  obscurants 
when  the  ratios  of  the  particle  size  to  the  wave¬ 
length  are  very  small  and  very  large,  respectively. 
When  the  wavelength  of  the  incident  radiation  ap¬ 
proaches  the  size  of  the  interacting  particles,  the 
mostly  widely  employed  solution  technique  is  Mie 
theory,  which  is  an  exact  solution  to  Maxwell’s 
electromagnetic  theory  for  a  plane  wave  incident 
on  a  homogeneous  sphere.  Kumai  and  Russell 
(1969)  and  Seagraves  (1981)  applied  Mie  theory  to 
determine  the  obscuring  effects  of  ice  fog  for 
various  wavelengths  in  the  visible  and  infrared 
regions.  Munis  and  Delaney  (1972)  and  O’Brien 
and  Kumai  (1973)  conducted  transmission  studies 
in  laboratory-generated  ice  fog  and  found  general 
agreement  between  the  experimental  extinction 
values  and  the  values  calculated  using  the  Mie 
theory. 


Application  of  Mie  theory 
to  ice  fog 

If  the  shape  of  the  ice  fog  particles  is  assumed  to 
be  spherical,  the  scattering  and  absorption  param¬ 
eters  of  ice  fog  in  the  visible  and  infrared  wave¬ 
lengths  can  be  derived  from  Mie  theory.  The  de¬ 
tailed  derivation  of  Mie  theory  can  be  found  in  the 
literature  and  is  not  repeated  here.  Only  the  com¬ 
monly  used  parameters  associated  with  Mie  theory 
are  defined  here.and  some  thoughts  on  the  applic¬ 
ability  of  Mie  theory  to  the  ice  fog  propagation 
model  are  presented. 

The  total  energy  scattered  by  a  panicle  is  derived 
by  integrating  the  scattered  energy  over  all  scatter¬ 
ing  angles.  The  scattering  cross-sectional  area  C„„ 
is  then  defined  as  the  area  that  intersects  as  much 
radiation  as  is  scattered  by  the  particle.  The  ex¬ 
pression  of  the  energy  scattered  in  terms  of  area 
permits  normalization  of  the  scattering  parameter. 
The  normalized  scattering  parameter  commonly 
used  is  the  scattering  efficiency  factor,  which  is  the 
ratio  of  the  scattering  cross-sectional  area  to  the 
geometrical  cross-sectional  area  of  the  scattering 
particle.  The  Mie  efficiency  factors  for  extinction 
and  absorption  are  normalized  in  the  same  man¬ 
ner.  The  efficiency  factors 

Q,ca  =  C.cJirr1,  Q.„,  =  C.t./rr1, 

and  Qr„  “  CMt/rr 

for  ice  fog  particles  in  the  visible  and  infrared 
regions  calculated  using  a  computer  program  for 
Mie  scattering,  AGAUS  82  (Miller  1983),  are 
shown  in  Figure  4.  The  complex  index  of  refrac¬ 
tion  of  ice  used  in  the  calculation  is  from  Irvine 
and  Pollack  (1968). 

The  scattering  cross-sectional  area  discussed 
previously  does  not  contain  information  regarding 
angular  scattering.  The  intensity  distribution  of 
scattered  energy  as  a  function  of  scattering  angle  is 
referred  to  as  the  phase  function.  This  function  is 
integrated  over  the  detector’s  acceptance  solid 
angle  in  order  to  correct  for  the  forward-scattered 
energy  that  enters  the  field  of  view  during 
transmittance  measurement. 

The  applicability  of  Mie  theory  for  use  in  an  ice 
fog  transmission  model  requires  experimental  data 
on  the  optical  properties  of  ice  fog  to  support  the 
equivalent-sphere  model  presumed  by  Mie  theory. 
Thuman  et  al.  (1955),  in  an  effort  to  determine  the 
physical  state  of  the  suspended  ice  fog  particles, 
measured  the  intensity  of  visible  light  scattered  by 
ice  fog  in  situ  on  10  different  occasions  at 
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Scattering  Angle 

Figure  5.  Comparison  of  experimental  and  calculated 
phase  functions. 

temperatures  ranging  from  -34°  to  -40°C.  The  more  closely  resembles  scattering  patterns  ob- 

scattered  intensity  was  measured  in  a  horizontal  served  for  irregular  particles.  An  important  obser- 

planeat  10°  intervals  from  15  to  165°.  Concurrent  vation  not  shown  in  Figure  5  was  that  no  signifi- 

observation  of  the  settled  ice  fog  particles  in-  cant  differences  were  apparent  in  the  light-scatter- 

dicated  that  between  85%  and  96%  of  the  par-  ing  properties  of  the  droxtal  types  and  the  more 

tides  were  of  the  droxtal  type.  (The  large  percen-  developed  ice  fog  crystals.  This  implies  that  the 

tage  of  droxtals  at  the  reported  temperature  range  same  phase  function  can  be  used  to  describe  the 

is  not  explained.)  This  large  number  of  spherical  scattering  properties  of  the  different  types  of  ice 

particles  provided  a  fortuitous  opportunity  for  fog  in  the  visible  region. 

comparing  the  measured  and  calculated  values  of  The  limited  data  available  for  the  radiation  scat- 

the  angular  variation  of  scattered  light.  An  addi-  tering  properties  of  ice  fog  necessitate  more  defini¬ 
tional  opportunity  for  comparative  analysis  was  live  experiments  in  the  future.  Until  these  meas- 

provided  by  Thuman  with  his  scattering  measure-  urements  are  available,  phase  function  measure¬ 
ment  conducted  for  ice  fog  composed  of  mostly  ments  conducted  for  artificially  generated  ice  crys- 

prismatic  columns  and  hexagonal  plates.  The  tals  can  be  evaluated  for  analogous  results.  Lab- 

■<?an  radii  reported  for  the  ice  fog  composed  of  oratory  measurements  of  the  phase  function  of 

ly  droxtals  and  the  more  developed  crystals  laboratory-generated  ice  particles  provide  a  realis- 

6.0  and  6.5  pm,  respectively.  tic  representation  of  an  in-situ  ice  fog.  These 

AGAUS  82  computer  program  for  Mie  measurements  (Huffman  and  Thursby  1970  and 

c  .ring  (Miller  1983)  was  utilized  to  obtain  the  Sasson  and  Liou  1978)  again  indicated  greater  side 

leoretical  scattering  pattern  for  an  ensemble  of  scattering  by  ice  crystals  than  indicated  by  the 

;e  spheres  with  similar  size  distribution  to  the  ice  computed  results  for  spherical  assemblies.  Given 

og  encountered  by  Thuman.  Comparison  of  the  the  experimental  evidence  accumulated  thus  far, 

measured  and  calculated  values  normalized  to  one  can  conclude  that  Mie  theory  can  be  applied 

coincide  at  40°  is  shown  in  Figure  5.  The  presence  for  forward-scattering  problems  in  ice  fog,  but 

of  a  minimum  and  a  maximum  at  the  respective  limitations  exist  for  side  and  backscatter  calcula- 

lateral  and  backscatter  angles  in  the  theoretical  tions. 

values  does  not  occur  with  the  experimental  re-  In-situ  measurements  of  infrared  radiation  scat- 

suits.  Rather,  the  shape  of  the  experimental  values  tering  by  ice  fog  do  not  appear  in  the  literature, 

for  the  angular  variation  of  scatter  radiation  However,  laboratory  measurements  of  infrared 

departs  from  the  equivalent  sphere  model  and  (10.6  pm)  scattering  by  artificially  generated  ice 


6 


Rearranging  this  equation  for  a. 


crystals  of  plate  and  columnar  structure  provide 
valuable  insight  into  the  potential  scattering  prop¬ 
erties  of  ice  fog  in  the  infrared  region  (S?sson 
1981).  The  angular  scattering  pattern  produced  by 
the  ice  crystals  again  deviates  from  Mie  values  em¬ 
ploying  the  equivalent-sphere  model.  As  in  the  vis¬ 
ible  regions  where  the  near-forward  scattering  ef¬ 
ficiency  can  be  determined  from  Mie  theory,  for¬ 
ward-scattering  by  ice  fog  in  the  infrared  regions 
may  also  be  determined  in  this  manner. 

The  experimental  data  available  appear  to  indi¬ 
cate  that  the  utility  of  Mie  theory  for  an  ice  fog 
propagation  model  in  the  visible  and  infrared  re¬ 
gions  of  the  spectrum  is  limited  to  forward  scatter¬ 
ing.  Applications  requiring  knowledge  of  the  scat¬ 
tering  properties  of  ice  fog  in  all  directions  will  re¬ 
quire  modeling  of  the  scattering  process  employ¬ 
ing  judicious  use  of  empirical  and  theoretical  re¬ 
sults. 


indicates  that  the  value  for  a  is  determined  by  N. 
The  extinction  coefficient  for  an  ice  fog  distribu¬ 
tion,  given  the  number  concentration,  can  now  be 
expressed  as 

n  -  N7  7  a 


? max 


Qrx,r  e\p[(-a/y)(r/rc)y)dr.  (7) 


Comparisons  of  infrared  to  visible 
volume  extinction  coefficients 

Comparison  of  the  transmission  properties  of 
various  wavelengths  through  an  obscuring  medi¬ 
um  can  serve  as  a  guideline  in  selecting  the  electro- 
optical  system  capable  of  optimal  performance 
through  the  given  medium.  With  this  in  mind,  the 
relationship  between  the  visible  and  infrared  ex¬ 
tinctions  through  ice  fog  is  investigated.  The  influ¬ 
ence  of  ice  fog  particle  size  on  extinction  proper¬ 
ties  in  the  visible  and  infrared  regions  has  been 
previously  illustrated  (Fig.  4).  Given  this  fact,  one 
can  conclude  that  the  shape  of  the  ice  fog  size  dis¬ 
tribution  curve  will  have  a  significant  role  in  influ¬ 
encing  the  infrared  to  visible  extinction  ratios.  In 
an  effort  to  determine  the  effect  of  the  size  distri¬ 
bution  in  this  relationship,  comparative  extinction 
values  were  calculated  for  a  broad  range  of  ice  fog 
distributions  expressed  in  terms  of  the  modified 
gamma  distribution. 

The  solution  for  determining  the  extinction  co¬ 
efficient  through  an  ice  fog  medium  represented 
by  the  modified  gamma  distribution  is  derived 
from  eq  1  and  3.  Forward-scattering  corrections 
(eq  4)  are  not  considered  in  the  following  calcula¬ 
tions  since  the  transmission  geometry  can  vary  de¬ 
pending  on  the  electro-optical  system  in  use.  Inte¬ 
gration  of  eq  1  (Dwight  1961)  over  the  entire  range 
of  radii  yields  the  total  number  of  particles  N  per 
unit  volume  where 


N  =  (a/y) 


a  , -(or  +  l)/> 
yr7  ) 


(5) 


The  volume  extinction  coefficients  were  calcu¬ 
lated  for  0.6328-Atm,  3.5-/un  and  10.6-pm  regions. 
The  size  distribution  parameters  were  chosen  to 
encompass  a  broad  range  of  distribution  shapes  in 
order  to  identify  the  influence  of  size  distribution 
on  the  extinction  properties  at  the  selected  wave¬ 
lengths.  The  relative  extinction  properties  were 
evaluated  for  115  different  distributions  repre¬ 
sented  by  the  following  combinations  of  parame¬ 
ters: 

or  =  2.0 

y  =  0.75,  1.0,  1,5,  1.95,  2.75 

rc  =  1.0-12.0  fim  at  0.5-/xm  intervals. 

These  combinations  of  parameters  do  not  neces¬ 
sarily  represent  actual  experimental  ice  fog  distri¬ 
bution  data.  They  are  presented  for  continuity  to 
demonstrate  the  sensitivity  of  the  infrared  to  visi¬ 
ble  extinction  ratios  to  the  modified  gamma  pa¬ 
rameters.  However,  the  distributions  for  ice  fog 
occurring  at  different  temperature  ranges  reported 
by  Ohtake  (1970)  have  been  included  in  this  inves¬ 
tigation.  One  can  also  assume  that  many  of  the 
distributions  are  in  fact  realistic  representations  of 
ice  fog  distributions. 

The  relative  performance  of  the  3.5-^m  wave¬ 
length  with  respect  to  the  visible  radiation  is  not 
favorable  (Fig.  6).  Although  the  transmission  at 
3.5  ^m  improves  with  increasing  particle  size  dis¬ 
tribution,  the  visible  radiation  is  less  obscured 
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Figure  6.  3 . 5 -pm  /0. 63 2 8- pm  extincticn  ratio. 


than  that  at  3.5  pm  for  all  distributions  likely  to  be 
encountered  in  ice  fog.  The  comparison  of  10.6  pm 
to  visible  radiation  is  perhaps  more  interesting, 
due  to  the  crossover  in  the  relative  performances 
of  the  two  wavelengths  (Fig.  7).  The  10.6-pm 
wavelength,  relative  to  the  visible,  degrades  no¬ 
ticeably  with  increasing  size  distributions,  indicat¬ 
ing  that  with  increasing  temperature  one  can  ex¬ 
pect  visible  systems  to  outperform  infrared  sys¬ 
tems.  Specifically,  the  extinction  ratio  of  the  10.6 
pm  to  0.6328  pm  in  ice  fog  formed  at  temperatures 
below  -38  °C  (e.g.  A,  representing  the  composite 
ice  fog  size  distribution  formed  around  -39°  to 
-41  °C  as  reported  by  Ohtake)  is  approximately 
0.55,  yet  actually  exceeds  unity  for  ice  fog  size  dis¬ 
tributions  associated  with  increasing  temperatures 
(e.g.  B,  showing  ice  fog  typical  of  that  formed 
around  -35  °  to  -36.9  °C). 

10.6-pm  extinction  as  a  function  of 
ice  fog  mass  concentration 

The  volume  extinction  coefficient  for  the  10.6- 
pm  wavelength  can,  within  limits,  be  approxi¬ 
mated  as  a  function  of  the  ice  fog  mass  concentra¬ 
tion  (A/C)  pseudo-independently  of  the  ice  fog  dis¬ 
tribution.  This  relationship  assumes  that  the  maxi¬ 
mum  radius  of  the  ice  fog  particles  in  any  given 
distribution  is  around  12.0  pm,  or,  if  larger  parti¬ 
cles  exist,  that  they  do  not  contribute  significantly 
to  either  the  B„,  or  the  A/C.  The  A/C  is  often  a 
more  easily  obtainable  parameter  than  the  size  dis¬ 
tribution.  Therefore  the  utility  of  this  relatively 
simple  relationship  can  be  incorporated  into  the 


B.  (a*2.y*l.95,re*8.5) 


(a  =  2,y  =2.75,rc*4.5) 
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Figure  7.  10. 6- pm/0. 6328-y.m  extinction  ratio. 


ice  fog  propagation  model.  The  procedure  for 
identifying  a  specific  wavelength  whose  extinction 
coefficient  is  a  function  of  A/C  was  outlined  by 
Chjflek  (1978)  for  fog  and  cloud  droplets. 

The  AfC  of  ice  fog  is  expressed  as 


'  max 

=  J  *r5  n(r)  dr 


where  e,(0.91  g/cmJ)  is  the  density  of  ice.  The 
volume  extinction  coefficient  is  related  to  the  A/C 
in  the  following  form: 


D  _  ^ 

-  A 


3  A/C  ^min 


'max 

J~  xr2  n(r)  Q.x,  dr 


'max 

j~  *rJ  n(r)  dr 


The  extinction  efficiency  factor  Q.„  for  the 
10.6-pm  wavelength  can  be  approximated  by  a 
linear  relationship  (Fig.  8),  such  that 

Q...  =  0.21  r  (10) 

where  the  maximum  radius  is  around  12.0  pm. 
Substituting  this  new  value  for  Q.„  into  eq  9  yields 
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Figure  8.  Linear  approximation  for  Q. 


B.„  =  (0.21) 


3  MC  r, 


'mu 

J~  rrl  n(r )  dr 


'max 

j~  *r3  n(r)  dr 


Finally,  canceling  of  the  integral  on  the  right  side 
of  the  equation  results  in  a  size-distribution-inde¬ 
pendent  relationship  between  the  10.6-/tm  extinc¬ 
tion  and  the  MC  of  ice  fog,  which  is  simply 

B.„  =  0.17  MC.  (12) 

In  order  to  evaluate  how  accurately  the  MC- de¬ 
pendent  extinction  values  approximate  the  exact 
solution,  comparative  results  were  analyzed.  The 
composite  size  distributions  for  the  three  tempera¬ 
ture  ranges  reported  by  Ohtake  (Fig.  2)  were  eval¬ 
uated.  The  MC  values  of  the  three  distributions 
for  the  arbitrarily  selected  number  density  of  100 
particles/cm  *  were  derived  by  multiplying  n(r)  in 
eq  1  by  %  rr’  q  ,  and  integrating  to  yield 
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Figure  9.  10.6-nm  extinction  as  a  function 
of  ice  fog  mass  concentration. 


Substituting  the  value  a  in  eq  6  gives 


Afr  ",  v/  «  V3/>  ri(a+4)/'yl 
A/C=/txe,N(— ) 


Then,  using  the  approximate  relationship  (eq  12) 
the  volume  extinction  coefficients  were  calculated. 
Also,  the  extinction  coefficients  were  calculated 
using  the  exact  solution  (eq  7)  for  identical  num¬ 
ber  density.  The  comparative  results  are  shown  in 
Figure  9.  The  line  represents  the  extinction  coeffi¬ 
cient  as  a  function  of  the  ice  fog  mass  concentra¬ 
tion  independent  of  the  size  distribution.  The 
result  seems  to  indicate  that,  given  only  the  MC  of 
ice  fog,  the  10.6-pm  extinction  coefficient  can  be 
approximated  within  reasonable  accuracy.  Addi¬ 
tional  verification  of  the  MC  to  10.6-/un  extinc¬ 
tion  relationship  is  demonstrated  by  similar  com¬ 
parative  analyses  of  size  and  mass  distributions  of 
ice  fog  (also  shown  in  Fig.  9)  reported  by  Kumai 
(1964). 


CONCLUSIONS 

Ice  fog  can  present  serious  visibility/electro- 
optical  system  performance  problems  during  mili¬ 
tary  operations  in  arctic  and  subarctic  regions. 


n;n.r,rer 


Thus,  knowledge  of  the  physical  and  optical  prop¬ 
erties  of  ice  fog  is  necessary  to  predict  the  perfor¬ 
mance  of  electro-optical  systems  in  this  environ¬ 
ment.  It  appears  that  the  potential  exists  not  only 
for  predicting  the  formation  of  ice  fog,  but  also 
for  determining  its  physical  properties,  given  the 
nature  of  the  water  vapor  source  and  the  ambient 
conditions. 

Based  on  experimental  evidence  to  date,  it  may 
be  concluded  that  Mie  theory  can  be  applied  to  de¬ 
termine  forward-scattering  properties  of  ice  fog  in 
the  visible  and  infrared  wavelengths,  but  that  lim¬ 
itations  exist  for  applications  involving  significant 
side  and  backscatter. 

Visible  radiation  is  somewhat  less  subject  to  ob¬ 
scuration  by  ice  fog  than  is  infrared  radiation  of 
3.5-pm  wavelength,  regardless  of  temperature. 
However,  the  extinction  of  10.6-pm  radiation  rela¬ 
tive  to  that  of  visible  radiation  is  dependent  on  the 
size  distribution  of  the  ice  fog,  which  is  in  turn 
temperature-related.  It  appears  that  ice  fog  is  con¬ 
siderably  more  transparent  to  10.6-pm  radiation 
than  to  visible  radiation  at  temperatures  below 
about  -38  °C,  whereas  at  higher  temperatures  the 
reverse  is  likely  to  be  true. 

On  the  basis  of  some  reasonable  assumptions 
concerning  the  probable  lack  of  contribution  to 
obscuration  by  large  ice  particles  in  ice  fog,  it  has 
been  shown  that  the  volume  extinction  coefficient 
at  10.6  /im  may  be  represented  as  being  directly 
proportional  to  the  mass  concentration  of  ice  fog 
particles  in  the  air.  The  proportionality  constant  is 
shown  to  be  0. 17  MC. 
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